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Many-Body Perturbation Theory: 
The GW Approximation 

Christoph Friedrich and Arno Schindlmayr 

Institute of Solid State Research 
Forschungszentrum Jiilich 
52425 Jiilich, Germany 
E-mail: {c.friedrich, a. schindlmayr} @fz- juelich.de 

In this lecture we present many-body perturbation theory as a method to determine quasiparticle 
excitations in solids, especially electronic band structures, accurately from first principles. The 
main ingredient is the electronic self-energy that, in principle, contains all many-body exchange 
and correlation effects beyond the Hartree potential. As its exact mathematical expression is 
unknown, approximations must be used in practical calculations. The GW approximation is 
obtained using a systematic algebraic approach on the basis of Green function techniques. It 
constitutes an expansion of the self-energy up to linear order in the screened Coulomb poten- 
tial, which describes the interaction between the quasiparticles and includes dynamic screening 
through the creation of exchange-correlation holes around the bare particles. The implemen- 
tation of the GW approximation relies on a perturbative treatment starting from density func- 
tional theory. Besides a detailed mathematical discussion we focus on the underlying physical 
concepts and show some illustrative applications. 

1 Introduction 

In the previous lectures we have seen that density functional theory (DFT) is the method 
of choice when we are interested in the ground-state properties of a many-electron sys- 
tem. DFT is based on the Hohenberg-Kohn theorem, 1 which states that there is (a) a 
one-to-one correspondence between the ground-state density no (r) and the external po- 
tential as well as (b) a variational principle for the energy functional E[no] < E[n]. The 
second statement allows to obtain the ground state of a many-electron system by variation 
of its density, a quantity that is much less complicated than the many-electron wave func- 
tion ^o( r i; • ■ ■ j r iv), where N is the particle number. The first statement implies that the 
many-particle Hamiltonian is a functional of the ground-state density. Since the diagonal- 
ization of the Hamiltonian yields the complete excitation spectrum, the excited states can 
ultimately be regarded as functionals of the ground-state density as well. However, the 
Hohenberg-Kohn theorem does not provide us with an explicit mathematical form. In this 
lecture we show that excited-state properties can be accessed more directly with a purpose- 
built method, the so-called many-body perturbation theory. 2 3 Incidentally, in practice its 
implementation within the GW approximation 4 for the electronic self-energy is based on a 
perturbative evaluation with Kohn-Sham orbitals and can, therefore, finally be interpreted 
as the desired density functional. 

The solution of the Kohn-Sham equation 5 of DFT yields a whole spectrum of single- 
particle states, and one is tempted to identify the corresponding eigenvalues with excitation 
energies. Strictly speaking, such an interpretation is wrong: the Kohn-Sham wave func- 
tions and eigenvalues must be considered as mathematical tools and cannot be endowed 
with a physical meaning. The only exception is the energy of the highest occupied state, 
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Figure 1. Compai'ison of LDA, GW and experimental band gaps for a vai'iety of materials. Taken from Ref. 



which equals the exact ionization potential (or chemical potential for metals). 6 ' 7 Con- 
sequently, while often qualitatively correct, the DFT band structure fails to give reliable 
quantitative values for the band gaps of insulators and semiconductors, which are often 
underestimated by as much as 1.0 eV or more. In the case of Ge the local-density approx- 
imation (LDA) of DFT even predicts a semi-metal with a negative band gap rather than 
a semiconductor. In this lecture we demonstrate that the Kohn-Sham eigenvalues can be 
corrected using Green function techniques and the GW approximation for the electronic 
self-energy. Figure 1 shows a comparison of LDA and self-energy corrected band gaps 
with respective experimental values for a variety of materials. The underestimation within 
the LDA as well as the improvement by the GW approximation are evident. 

Band gaps are experimentally measured by photoelectron spectroscopy. Figure 2 gives 
a schematic illustration. In direct photoelectron spectroscopy a photon with energy fou 
impinges on the sample and ejects an electron, whose kinetic energy £kin is subsequently 
measured. The binding energy of this electron is given by the difference = i?kin — hoj. 
Actually, we already simplified the argumentation here, as the formulation "binding energy 
of an electron" suggests that the electrons are independent. In reality they are correlated 
through the Coulomb interaction, and the ejection of an electron is always a many-body 
process. In this general sense equals the difference ej = Eq — E^ 1 between the 
total energy Eq of the TV-particle ground state and the energy E*- 1 of the (TV - 1)- 
particle state ^f^ 1 that remains after the emission. Inverse photoelectron spectroscopy is 
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Figure 2. Schematic illustration of direct and inverse photoelectron spectroscopy. In both processes the particle 
number changes. The measured energy difference iS^in — &<J corresponds to ej 
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the complementary process: electrons are injected into the sample, and the energy of the 
emitted photon is measured. The number of electrons in the system thus increases from TV 
to iV + 1, and we can identify £kin — fuo with the energy difference ej = E^ +1 — E$ of 
the many-electron systems. 

The fact that the independent-electron picture breaks down due to the strong Coulomb 
interaction questions single-electron concepts like band structure or Fermi surface. Still, in 
practice these work surprisingly well. In fact, we can at least retain a nearly-independent- 
particle picture if we consider quasiparticles instead of electrons (or holes). In the case 
of electron injection into a sample the repulsive Coulomb interaction creates a Coulomb 
hole around the additional electron (see Figure 3). Analogously, if an electron leaves the 
system, its Coulomb hole also disappears. Relative to the ground-state iV -electron system, 
the addition (removal) of an electron in indirect (direct) photoelectron spectroscopy hence 
creates (annihilates) an ensemble consisting of the bare electron and its oppositely charged 
Coulomb hole. This ensemble behaves in many ways like a single-particle and is thus 
called "quasiparticle". Since the Coulomb hole reduces the total charge of the quasipar- 
ticle, the effective interaction between quasiparticles is screened and considerably weaker 
than the bare Coulomb interaction between electrons. In fact, the screened interaction is 
sufficiently small so that the quasiparticles can be regarded as approximately independent, 
which finally justifies the independent-particle approximation and explains the success of 
mean-field theories. 

A theoretical description of processes involving the ejection or injection of electrons 
requires a framework that links the -particle with the (N ± l)-particle systems. For 
this purpose we employ many-body perturbation theory. The central variable is the time- 
ordered Green function G(rt, r't'). As we will see, it contains the excitation energies ej 
and even the excitation lifetimes. Besides, we can directly obtain the ground-state electron 
density, the expectation values of one-particle operators and the ground-state total energy 
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Figure 3. The electrons in a many-electron system are correlated by the strong Coulomb interaction v. The 
motion of one electron depends on the motion of all other electrons. A nearly-independent-particle picture can 
be recovered within the quasiparticle concept. Due to exchange and correlation effects a Coulomb hole forms 
around an electron and behaves together with it like a single entity, which is called quasiparticle. Quasiparticles 
interact via a weak screened interaction W instead of the strong Coulomb interaction. 

from it. The Green function is hence capable of giving access to the same observables as 
the ground-state electron density. In contrast to the DFT expression E [n], the functional 
E [G] is even known exactly. 2 While the Green function contains much more information 
than the electron density, it is also a more complicated function and thus rarely applied 
to ground-state properties. In the present lecture we will, therefore, concentrate on the 
calculation of excited states. 

Section 2 lays the theoretical foundations of the method. More complicated derivations 
are deferred to the appendix. The GW approximation is discussed in Section 3.1, and 
some aspects of its numerical implementation are given in Section 3.2. As an illustration, 
Section 3.3 presents a number of selected applications. Section 4 contains the summary. 

2 Theory 

2.1 Green Function 

In this section we introduce the time-ordered Green function and examine its properties. 
We use the second-quantization formulation of quantum mechanics. 2 3 For the present 
purpose it is sufficient to know that this formalism involves field operators ^(r) and V>^(r) 
that describe the annihilation and the creation of an electron at the position r, respectively. 
We will not take spin dependence explicitly into account. If necessary, the spin quantum 
number can simply be added to the formulas by considering it to be part of the spatial 
coordinate r. 

The Green function G e (rt,r't') is defined such that ihG e (rt,r't') is the probability 
amplitude for the propagation of an additional electron from (r', t 1 ) to (r, t) in a many- 
electron system with the Hamiltonian (Eq. 41)). This process brings the system from the 
TV-electron ground state \^(t')) to a final state tp(r)U(t, t')fo(r') (*'))• The final 
state is constructed by the successive action of the electron creation operator ^)t(r'), the 
evolution operator U(t, t') = exp[— iH(t — t')/H], which takes the system from t' to a later 
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time t > t', and the electron annihilation operator ip(r) on the TV-electron ground state. As 
the probability amplitude is given by the overlap of the final state with (£)), the Green 
function becomes 



where 6{t — t') is the Heaviside step function defined by 

^-*o ={;;; 



*«{*)) 0(t-t) 



ift > t' , 
iff < t' . 



(1) 



(2) 



For the last equality in Eq. (1) we changed from the Schrodinger to the Heisenberg picture, 
where the expression is particularly simple. States and operators in the two pictures are 
related by 

|* H >=#(0,f)|*s(f)> and An(t) = U(0,t)A a U(t,0). (3) 
In the following we always omit the indices S and H. Similarly, we have the Green function 

G h (T't',Tt) = - l - \ft(T't')i>(Tt)\v»)6(t' - t) (4) 



for the propagation of an additional hole from (r, t) to (r', t') . As a matter of convenience, 
we combine the two expressions in one time-ordered Green function 

G{vt,v't') = G e (rt,r't') - G h (r't',rt) = - % - \f [^(rf^VO] | *o) , (5) 

where we used the time-ordering operator T, which rearranges a series of field operators 
in order of ascending time arguments from right to left with a factor (—1) for each pair 
permutation. Depending on the time order, Eq. (5) describes either electron (t > t') or 
hole (t < t') propagation. The electron density n(r) can be expressed in terms of the 
Green function as 



n(rt) = ft(rt)$(rt) = -ihG(rt,rt + rj) 



(6) 



Here and in the following r) is an infinitesimal positive number. It serves only to enforce 
the correct order of the field operators. Its unit should always be clear from the context; 
presently it is an infinitesimal time. 

Let us consider the time-ordered Green function G(r, r'; r) of a stationary system with 
r = t — t'. If we insert the closure relation £\ l^^ 1 ) (\E r ^ v±1 1 = 1 between the two 
field operators in Eq. (5), where {l^l/f^ 1 )} is the complete set of state vectors of the 
(N ± l)-particle system, transform to the Schrodinger picture and use the definitions 



Vf- 1 (r) = (*f- 1 |vi(r)|<) and Vf +1 (r) = |^(r) 



(7) 



together with the excitation energies 



JV-l eiJV jpN-l 



and 



-E, 



N 
J 



(8) 



5 



then we obtain 



G(r,r';r) = ^ W +1 V)e" fc " 



+ ^E^ 1 (r)V-r i *(r') e - e "" V ^(-r) . (9) 

i 

The sums run over the ground state and all excited states of the (N — 1)- and (N + 1)- 
particle system, respectively. Expression (Eq. (9)) can be interpreted as follows: The state 
after the addition of an electron (r > 0) is represented by a linear combination of excited 
states 

^ (r ')|^) = ^^+ 1 *(r')|*f +1 ) (10) 

i 

that subsequently evolve according to their respective phase factors exp(— ie? +1 t/H). The 
resulting state is then probed at the point r by the projections ipf +1 (r). The case r < 
(hole propagation) is analogous. Consequently, the Green function indeed contains the 
complete excitation spectrum of the (TV ± 1) -particle system. Fourier transformation of 
Eq. (9) to the frequency axis using the Fourier transform of the Heaviside step function 

1 f°° i 
0( W ) = — / 0(t) e iwT -^dr = - . — , (1 1) 

finally yields the Lehmann representation of the Green function 

^.^)-E r (r C (r) + ^ i i(r ^" (r) - <i2 » 

^-r 1 hu - e i ^ + ir) ^ hui-e i -it] 

We observe that the Green function has poles at the true many-particle excitation energies 
£ jv±i T nese energies correspond to excitations of an (N — l)-particle and an (N + 1)- 
particle system and hence to those processes measured in direct and inverse photoelectron 
spectroscopy. In the case of a noninteracting (or mean-field) system the ■0^ v+1 (r) are 
simply the unoccupied and the ^>^ v_1 (r) the occupied single-particle wave functions, the 
e^ 1 are the corresponding single-particle energies. In order not to overload the notation, 
we will drop the (N ± 1) superscripts from now on. 



2.2 Spectral Function 

In connection with Eq. (9) we can define the spectral function A(r, r';ui), i.e., the density 
of the excited (or quasiparticle) states that contribute to the electron or hole propagation. 
In a finite system this density is a series of delta functions at the excitation energies 

A(t,t';u) = ^ i (r)^(r')<S(foj-e i ) (13) 

i 

weighted by the products of the corresponding projections (Eq. (7)). This allows us to 
rewrite Eq. (12) as an integral over frequencies 

^(r, r,w) ^ J ^ - fy^i -|- S gn(fej' - n) iif^ 



6 



finite system 



A(m) 
A 



e,€lR 



position 



infinite system 




quasiparticle position -width 

— > lifetime 



Figure 4. The excitation peaks of a finite system in the spectral function A(u>) merge into quasiparticle peaks 
of finite width in the case of an infinite system. This gives rise to finite excitation lifetimes determined by the 
inverse of the peak widths. 



with max (e^ ) < fi < min (e£ +1 ) . In an infinite system corresponds to the chemical 
potential. The inequality max (e^ -1 ) < min (e^ +1 ) follows from the convexity of the 
total energy as a function of the particle number, i.e., Eq -1 — E$ > E$ — E$ +l : we lose 
more energy when removing an electron than we gain by adding one. With the identity 



— L_ = V ( - ] ± ind(x) 

X^IT] \X ' 



in the limit rj — > + , where V(l/x) is the principal value of l/x, we find that 
A(t, r'; ui) = —sgn(hu — fx) — Im G(r, r'; ui) . 

7T 

The closure relation of the functions in Eq. (7) yields another important property 

/oo 
A(r, r'-u)dw = £ ^(r)V* (r') = d(r - r') . 
-OO 



(15) 



(16) 



(17) 



When we change from a finite to an infinite system, the delta functions in A(r, r'; u>) merge 
and form a series of smooth peaks with finite line widths instead of sharp resonances (see 
Figure 4). However, if the resulting spectral features are of Lorentzian form, i.e., 

A(r,r';u>) = £>(r)Vn r ')- ^V^7 , (18) 



where the are the peak positions and |Tj| the corresponding peak widths, then we can 
perform the integration in Eq. (14) analytically and again obtain a discrete sum over i as in 
Eq. (12), provided that the energies are defined as complex numbers = + iTj. Conse- 
quently, the form of the Fourier transform (Eq. (9)) remains unchanged, too. The imaginary 
component of leads to a damping term exp(— \TiT \ / h), revealing that the excitation has 
a finite lifetime of ^ | . Physically, the de-excitation proceeds via Auger transitions 
that create electron-hole pairs on the way. The damping of the particle propagator G may 
seem surprising, as it suggests that the particle gradually disappears. However, one must 
keep in mind that we deal with an infinite system, i.e., N — > oo, and an additional particle 
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(electron or hole) can "dissipate" into the Fermi sea. In this sense, one often speaks of 
finite quasiparticle lifetimes and calls Vi( r ) an d e i the quasiparticle wave functions and 
energies, respectively. The quasiparticle equation (22) introduced in the next section holds 
for infinite systems if one uses an analytic continuation of the self-energy into the complex 
frequency plane. 

2.3 Dyson Equation 

Appendix A shows that the time-ordered Green function G(r, t';cj) of the interacting sys- 
tem obeys an integral equation, the Dyson equation 

G(r,r';w) = G (r, r'; w) + J J Go(r,r»E(r",r'»G(r"',r'; W ) ( i 3 r" ( i 3 r"', (19) 

where Go (r, r'; w) is the Green function of a mean-field system defined by 

fto¥>?(r)=eM(r) (20) 
with the single-particle Hamiltonian 

fto(r) = -|^V 2 + y ext (r) + "7~~~ f j^dV . (21) 
2m 47T£ J \r — r'| 

The quantities V e xt( r )> m > e an d £ o are defined as in Eq. (41). The Green function 
Go(r, r'; ui) is obtained from Eq. (12) with the wave functions ip® (r) and energies e?. The 
nonlocal and frequency-dependent function S(r, r',w) is the non-Hermitian self-energy 
operator, which contains all many-body exchange and correlation effects beyond the elec- 
trostatic Hartree potential. This can be more easily seen in a reformulation of the Dyson 
equation. By inserting the Lehmann representation (Eq. (12)) into Eq. (19), we find that 
the wave functions tpi (r) and energies obey the quasiparticle equation 

fto(rM(r) + J S ( r > r '; e*//i)V>i(r')dV = e^i(r) (22) 

(see Appendix B), which is nonlinear in e^. Although it looks very similar to the one- 
particle equations of mean-field approaches like Hartree, Hartree-Fock or DFT, it does not 
constitute a mean-field formulation, since the self-energy takes all dynamic many-electron 
processes into account. Consequently, the functions ^j(r) and energies must not be 
understood as single-particle quantities. In fact, they are defined in Eqs. (7) and (8) as 
properties of the many-electron system. From the nonlinearity of the quasiparticle equation 
it follows that the wave functions tpi (r) are not orthonormal, in contrast to single-particle 
wave functions. However, they do fulfill the closure relation (Eq. (17)). 

The Dyson equation (19) can be rewritten in the form of a geometric series by subse- 
quently replacing G on the right-hand side by Go + GoSG, which leads to, symbolically 
written, 

G = Go + GoSGo + GoSGoSGo + GoSGoSGoSGo + . . . . (23) 

This is a typical equation of scattering theory, where the different terms of the geometric 
series describe single, double, triple, etc., scattering processes, and S is the scattering po- 
tential. Such a succession of scattering processes can be illustrated by Feynman diagrams, 



8 




r 



Figure 5. Illustration of a series of scattering processes using Feynman diagrams. All zigzag lines representing 
the instantaneous Coulomb interaction must be drawn horizontally. Arrows going forward in time represent 
electron and those going backward in time hole propagators. The self-energy is the sum of all possible single 
scattering processes. 



where Go is drawn as a straight arrow and the Coulomb interaction as a zigzag line. Ac- 
cording to Eq. (23), a diagrammatic representation of a multiple scattering process should 
involve a series of arrows (Go) divided by single scattering processes (£). In the example 
of Figure 5 these are the exchange interaction, the creation of an electron-hole pair (the 
"bubble" diagram) and finally the creation of a pair that itself creates another pair. In order 
to obtain the complete Green function, we have to sum all multiple scattering processes, 
of which the one shown in Figure 5 is merely one example. The self-energy is given by 
the sum of all single scattering processes. The interpretation in terms of scattering pro- 
cesses allows to construct approximations for £ by the summation of diagrams considered 
essential for the physical behavior of a given electron system. In general, however, such 
approximations are rarely convergent, and too many processes turn out to be quantitatively 
important. Therefore, we apply a systematic algebraic method instead. 



3 Implementation and Applications 
3.1 GW Approximation 

In practice we must use an approximation for the self-energy, such as the GW approxima- 
tion, which contains the electron exchange and a large part of the electron correlation. It 
is formally derived in Appendix A and has a very simple mathematical form in the time 
domain 

E GW (r, r'; r) = ihG (r, r'; r)W(r, r'; r + 7?) . (24) 
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In order to calculate the self-energy contribution to the quasiparticle energies, we need the 
Fourier transform on the frequency axis 

S Glv (r,r';w) = ^ / G (r, r'; w + cj')W(t, r'; cj')e iu W . (25) 

2 7i" J-oo 

The first function on the right-hand side is the Green function of the noninteracting sys- 
tem defined by Eq. (20) and the second function the dynamically screened interaction 
W(r, r';ui), which is related to the bare Coulomb potential v(r, r') = e 2 / (4iveo |r — r'|) 
through the inverse of the dielectric function 

W(r, r'; u) = J e' 1 (r, r"; u)v{r", r')d 3 r" = v(r, r') + J n ind (r, r"; u)v(r", r')d 3 r" . 

(26) 

The screened interaction W(r, r'; w) is the effective potential at r' induced by a quasipar- 
ticle at r: the Coulomb potential of the electron repels other electrons in its neighborhood 
and thus gives rise to the formation of an exchange and correlation hole, whose effective 
positive charge ni n <j(r, r"; w) screens the bare Coulomb potential v(r, r') (see Figure 6). 
Analogously, an effective negative charge screens the Coulomb potential of a hole. The 
screened interaction is considerably weaker than the bare Coulomb interaction. The GW 
approximation uses the random-phase approximation (RPA) 

e(r, r'; w) = 5(r - r') - J v(r, r")P(r", r'; u)d 3 r" , (27) 

P(t,t';t) = -i/iGo(r,r';r)Go(r',r;-r). (28) 

It corresponds to a subset of scattering processes in the many-electron system. Some of the 
respective diagrams are just the ones shown in Figure 5. Using expression in Eq. (12) for 
the Green function Go of the noninteracting system we observe that the Fourier transform 
of the polarization function in Eq. (28) is given by 

occ. unocc. 

P(r,r'; W )=]T ]T ^ (r) (r) cpf (r') ^ (r') 

i 3 

x I 1 l _ \ (29) 

yhu +e? - +irj hu - + - if) J 
in terms of the wave functions tp® (r) and energies . 




Figure 6. The formation of the Coulomb hole around an electron at r screens its Coulomb potential v(r, r'). 
This leads to the definition of the screened interaction W(r, r') that takes into account the combined potentials 
of the bare electron and its screening cloud 7ii n d. The ensemble consisting of the electron and its polarization 
cloud is called "quasiparticle". 
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The well-known Hartree-Fock equations can be recovered from Eq. (22) if we use the 
energy-independent self-energy 

S HF (r, r') = ihG (v, r'; -r))v(v, r') (30) 

(given in the time domain) instead. By comparison with Eq. (24), we see that the GW 
approximation constitutes an expansion of the self-energy up to first order in the screened 
interaction as opposed to the bare Coulomb interaction in Eq. (30). This approximates 
the exact self-energy considerably better, because W is much smaller than v. Due to 
the similarity of the two self-energy expressions, the GW approximation can formally be 
regarded as a Hartree-Fock approach with a dynamically screened interaction W instead 
of the static Coulomb interaction v. 



3.2 Numerical Implementation 

For band structure calculations it is more efficient to obtain the directly from the quasi- 
particle equation (22) instead of solving the Dyson integral equation (19) and searching 
for the poles of the Green function. Furthermore, it is then possible to exploit the formal 
similarity to the Kohn-Sham equation 

h <pf s (r) + Vxc(i>f S (r) = ef s <pf s (r) , (31) 

where V^ c (r) is the local exchange-correlation potential. In many cases the Kohn-Sham 
eigenvalues ef s already provide a reasonable estimate of the band structure and are in 
qualitative agreement with experiment. For systems where the quasiparticle wave functions 
are known, one also finds ipf s (r) « ipi (r) . 9 This observation indicates that the self-energy 
correction £ (r, r';ei/h) — V xc (r)5(r — r') is small and justifies the use of first-order 
perturbation theory to obtain approximate energies 

d a ef s + (iff 8 |S (ei/h) - Vxckf S ) • (32) 

A solution of this nonlinear equation still requires the knowledge of the frequency depen- 
dence of the self-energy, which is not known in general. Therefore, we use the linear 
expansion 

W ; /hi ~ W '.JCS/fcA i 6i ~ ^ d^K 1 *' r '; e f S /^) /ton 

Mr,r ;ei/n) w L(r,r ;e j /H) H ^ , (33) 

which leads to 

e, « ef s + Z t (^f S |S (ef s /h) - V xc \ ^f s ) . (34) 
The quasiparticle renormalization factor is given by 

1 (ef s /h) 



Zi= i- (vP 



and equals the quasiparticle weight 




(35) 



Zi = J \iPi{v)\ 2 d 3 r < 1. (36) 
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With the decomposition of W into the bare Coulomb interaction v and the remainder W — 
v, the GW self-energy (Eq. (24)) splits into exchange and correlation parts, symbolically 
written as 

Y, aw = ihG$ s W = ihG$ s v + ihG$ s (W -v) = + . (37) 

Instead of Go we use the Kohn-Sham Green function Gq S . After inserting this decompo- 
sition into Eq. (25), we must evaluate the convolutions 

S° w (r,r';w) = ^ /°° GjP(r, r'; u + w>(r, rV w 'W , (38a) 

27T J-oo 

S° w (r,r';w) = ^ / G? s (r,r';u+w')[W(r I i';w')-»(r I r')](4 1 /. (38b) 

The integral (Eq. (38a)) can be evaluated analytically and leads to the well-known expres- 
sion for the Hartree-Fock exchange term 

2 oca KS* ( } KS , ) KS* f n KS f n 

x * 1 1 ' ' 47r£o4^y |r — r'l 

(39) 

In general, the second convolution (Eq. (38b)) must be computed numerically. For this pur- 
pose the integration contour is usually deformed to the complex plane, where the analytical 
continuations of Go and W are smoother. 

Let the Kohn-Sham wave functions be represented in a basis {C«(r)}. According to 
Eq. (29) we can then write the polarization function and all related quantities in terms of 
products x M (r) = Ca( r )C^( r ) w i tn tne composite index /z = (a, 0) as 

P(t,i>;u) = ^P F ( W ) X ;(r)^(r') . (40) 

The Eqs. (24) to (28) are solved by matrix operations: 

1. A self-consistent DFT loop produces the Kohn-Sham wave functions ipf s (r) and 
energies ef s . At this point we can already evaluate the exchange term (Eq. (39)). 

2. The polarization matrix P p „(w) is calculated according to Eq. (29). 

3. The dielectric matrix is obtained from e M „(o;) = 5^ — ^ t; M7 P 7 y(w) and inverted. 

4. Next the screened interaction W M „(u;) = {w)v lv is calculated from a matrix 
multiplication of the inverse dielectric function with the Coulomb matrix. 

5. The correlation term (y?f s |E^ W '| ¥?f S ) is evaluated according to Eq. (38b) with a 
numerical contour integration on the complex frequency plane. 

6. Finally, approximate quasiparticle energies are obtained from Eqs. (34) and (35). 

The computation of the dielectric function, its inversion and the convolution (Eq. (38b)) 
are very time-consuming. Therefore, some (especially older) codes approximate the in- 
verse dielectric function by a so-called plasmon-pole model. 10 ' 11 These models replace the 
imaginary component of £ _1 (w), which has a peaked structure, by a sum of delta functions 
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at the corresponding frequencies. This simplification reduces the third step to a single ma- 
trix inversion of the static dielectric function at uj = and makes an analytic evaluation of 
the frequency integral (Eq. (38b)) possible. 

3.3 Examples 

Although Hedin's seminal article 4 was already published in 1965, it was not before the 
middle of the 1980s that the first ab initio calculations for real materials were reported in 
the literature. In spite of several approximations in the numerical treatment, which were 
necessary because of the lack of computer power back then, initial results were already 
very promising. Hybertsen and Louie 12 as well as Godby et al. 13 showed that the cal- 
culated band gap of Si fell within a margin of about 0.1 eV from the experimental value. 
Shortly afterwards the same authors reported band gaps for several other semiconducting 
materials that turned out to be equally accurate. 14 ' 15 After these pioneering studies the 
GW approximation was applied to a variety of semiconductors with great success (see, 
e.g., Figure 1). The principal effect of the GW self-energy correction on the band struc- 
ture of a semiconductor is to rigidly shift the valence bands up and the conduction bands 
down, thus opening the band gap. Figure 7 shows this effect for Si as an example. 

Not only the band gaps of semiconductors and insulators are improved by the GW self- 
energy correction, but the correlation-induced band narrowing of metals is also correctly 
described. The band narrowing reflects the higher effective mass of quasiparticles (the 
polarization cloud adds to the electron mass) compared to bare electrons. For this reason, 
the self-energy is sometimes also referred to as "mass operator". Figure 8 shows the energy 




l r x 



Figure 7. LDA band structure (dashed lines) of silicon with GW self-energy corrected valence and conduction 
bands (solid lines). The GW approximation shifts the corresponding bands up and down, respectively, but leaves 
the dispersion essentially unaffected. 
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Figure 8. Comparison of LDA (dashed), quasiparticle (solid line) and experimental (crosses) bands for Na. Taken 
from Ref. 16. 
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Figure 9. De-excitation dynamics measured in time-resolved two-photon photoemission spectroscopy (dia- 
monds) and calculated with the GW approximation (solid line). Taken from Ref. 17. 

dispersion of Na as an example. 16 The band narrowing brought about by the GW self- 
energy correction leads to nearly perfect agreement with experiment. 

The calculated excitation or quasiparticle lifetimes can be directly compared with two- 
photon photoemission spectroscopy. This experimental method allows to measure dynam- 
ical de-excitation processes in electronic systems. After a first photon has excited the 
electron system (creating a "hot" electron), a second photon probes the quasiparticle den- 
sity of states like in ordinary direct photoelectron spectroscopy. The time delay between 
the two photons can be tuned such that the system can be observed in different stages of the 
electronic de-excitation process. From a series of measurements one can thus deduce the 
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lifetime, which depends on the excitation energy, i.e., the energy of the first photon. In the 
example of Ag in Figure 9 the theoretical curve r = h \T\^ obtained from the imaginary 
parts of the quasiparticle energies (see Section 2.2) closely follows the experimental data 
points. 17 

4 Summary 

In this lecture we presented the GW approximation for the electronic self-energy, which 
allows to calculate excited-state properties like excitation energies and lifetimes. The 
self-energy describes scattering processes between electrons and, in principle, contains 
all exchange and correlation effects beyond the electrostatic Hartree potential. The GW 
approximation includes a subset of these scattering processes. Apart from exchange it de- 
scribes the creation of electron-hole pairs within the random-phase approximation (RPA) 
that leads to the formation of polarization clouds around the bare particles. The ensemble 
of an electron or a hole together with its polarization cloud behaves essentially like a sin- 
gle entity and is called a quasiparticle. The quasiparticles interact via a screened potential 
that is considerably weaker than the bare Coulomb interaction. This makes a perturbative 
treatment possible. In this respect, the GW approximation constitutes an expansion of the 
self-energy up to linear order in the screened interaction. It works well in a large class 
of systems where the polarization effects covered by the RPA play the dominant role in 
electron correlation, such as simple metals and semiconductors. 

The GW approximation is by nature a perturbative approach. Actual GW calculations 
are usually based on the self-consistent Kohn-Sham wave functions and energies as a start- 
ing point. This method has its limitations in materials where DFT already gives unphysical 
results. It breaks down for systems with very strong electronic correlation, which is in- 
sufficiently described by the available exchange-correlation functionals. The large error in 
the band gap of NiO in Figure 1 is an example. In reality, NiO is a strongly correlated 
Mott-Hubbard insulator, whereas it comes out as a semiconductor with a very small band 
gap (nearly a semi-metal) in DFT calculations. 

The GW method is designed for the analysis of excited states of the (N ± 1) -electron 
systems. The treatment of optical absorption processes, where the particle number does 
not change due to the promotion of valence electrons into unoccupied conduction states 
rather than emission, requires the simultaneous description of two particles, an electron 
and a hole, i.e., an exciton. Consequently, one must describe such a process with a two- 
particle Green function. In this case many-body perturbation theory leads to the so-called 
Bethe-Salpeter equation. Absorption spectra obtained from this equation are indeed very 
accurate. 18 An alternative is time-dependent density functional theory, 19 which also gives 
access to the excited states of an A^-electron system. 
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Appendix A 



Hedin equations 



With the field-operators introduced in Section 2. 1 we can rewrite the many-particle Hamil- 
tonian 



(41) 



where V^. x t( r ) is the potential created by the atomic nuclei, v(r, r') = e 2 / (47t£o |r — r'|) 
the Coulomb interaction, m the electron mass, e the electron charge and £ the vacuum 
dielectric constant, as 

H = J ft(r)h(r)$(r)d 3 r + \ j j ^ t (r)^ t (r')«(r,r')^(r')^(r)d 3 rd 3 r' (42) 

with the one-particle operator 

h 2 



Mr) = -^V 2 + Vext(r). 



(43) 



The expression in Eq. (42) is just a mathematical reformulation of Eq. (41) and should not 
be mistaken for the energy expectation value in Hartree theory, although it looks similar. 
From the equation of motion for the annihilation operator 

ih^(r,t) = [i>(r,t),H]_ =h(r)i>(r,t) + J v(r, r')^(r'> t)i>(r, t)d 3 r' , 

(44) 

which describes the time evolution of a Heisenberg operator in the same way as the 
Schrodinger equation describes that of a wave function, we can directly deduce the equa- 
tion of motion for the Green function 



ih^G(rt, r't') = S(r - r')S(t - t') + h(r)G(rt, r't') 



(45) 



i 



J „(r,r") \T [^(r",^(r",^(r.^ t (r'.0] | *") d 3 r" . 

This is not a closed equation, because it involves the two-particle Green function 

G(1234) = ~ (< \f [^(1)^(2)^(4)^(3)] | yty . (46) 

Here and in the following we denote the set of space-time coordinates (ri , t\) with a natural 
number 1, etc., and further define 



5(12) = <S(n - r 2 )6(h - h) 
v(12)=v(r 1 ,T 2 )6(t 1 -t 2 ), 

J dl = j d 3 n J dh , 
1 + = (ri.ti +v), 



(47) 
(48) 

(49) 
(50) 
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where T] is an infinitesimal positive time. With the two-particle Green function (Eq. (46)) 
we can rewrite Eq. (45) as 

ih^-G{\2) = 6(12) + ft(l)G(12) - ih / v (1+3) G (1323+) d3 . (51) 
oh J 

The additional infinitesimals in 1 + and 3 + make sure that the time order is the same as in 
Eq. (45). 

In order to employ the functional-derivative method, we introduce an external potential 
U(l) that is again set to zero at the end. Of course, all quantities from now on depend on 
U(l), while the equations remain invariant provided that we replace h(l) — > h(l) + U(l). 
We can use functional differentiation to define a number of useful quantities. The reac- 
tion of the density to changes in the external potential is governed by the linear-response 
function 

5n(l) 



R(12) 



(52) 

(7=0 



6U(2) 

The test potential and the Coulomb potential created by the induced charge can be com- 
bined into an effective potential 

U eS (1) = U(l) + jj v(13)R(32)U(2)d2 d3 , (53) 

which is related to U(l) via the inverse dielectric function 

= 5(12) + f v(13)R(32)d3 . (54) 

(7=0 J 



s -i (12)= «MD 



5U(2) 

With the definition of the polarization function 



P (12)= J%| (55) 



(7=0 



+ ... , 

(56) 



dU eS (2) 

and the chain rule for functional derivatives one obtains the geometric series 
£^(12) = 5(12) + J v(13)P(32)d3 + JJJ v(13)P(34)v(45)P(52)d3d4d5 

which can easily be inverted to yield 

£ (12) = 5(12) - J v(13)P(32)d3 . (57) 

If we take the Coulomb potential of an electron at 2 as the test potential, we get the screened 
potential 

W(12) = J £" 1 (13)^(32)d3 = v(12) + J J t;(13)P(34)W(42)d3d4 (58) 

as the effective potential at position 1. 

After this interlude we can go on with the derivation. For the functional derivative of 
the Green function with respect to the test potential we find 4 20 

5G(12) 



SU(3) 



= G(12)G(33+) - G(1323+) . (59) 

(7=0 
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This allows us to eliminate the two-particle Green function, and the integral in Eq. (51) 
becomes 

-ih J t;(l + 3)G(1323+)<i3 

= -ih (^J v{l3)G(33+)d3 S j G(12)+ih J ^(1+3) ^^ 3 (60) 

= V H (1)G(12) + J E(13)G(32)d3, 
where V H (1) is the Hartree potential [cf. Eq. (21)] and 

S(12) =ihjj v (l+3) 6 -^G- 1 (42)d3d4 

= -mJJ v(l + 3)G(14)^^-d3d4 (61) 

= ih JJ W(l + 3)G(14)T(42;3)d3d4 

the self-energy. For the second identity we used partial integration and for the third the 
chain rule for functional derivatives, the definition of the screened interaction (Eq. (58)) 
and the vertex function 

r(12;3) = 



SU eS (3) 



(62) 

(7=0 

With the self-energy (Eq. (61)) the equation of motion for the Green function (Eq. (51)) 
now becomes 



G(12) - J £(13)G(32)d3 = £(12) , (63) 

where we incorporated the Hartree potential into the one-particle operator 

ho(l)=h(l)+V H (l). (64) 

The delta function on the right-hand side of Eq. (63) demonstrates that G(12) is indeed 
a Green function in the mathematical sense. In a noninteracting system the self-energy 
vanishes, and Eq. (63) becomes 



i^-Ao(l) 



Go (12) =5(12). (65) 



Multiplication of Eq. (63) with Go and Eq. (65) with G from the left followed by integra- 
tion yields the Dyson equation 

G(12) = G (12) + JJ G (13)X(34)G(42)d3<i4. (66) 

Finally, Eqs. (65) and (66) allow us to rewrite the vertex function (Eq. (62)) as 

r(12; 3) = 5(12)5(13) + Jg-^, (67) 
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and with the identity 



S = ^// G(14)c " (45)0<52) ' i4 ' i5 

_ 2 5G{12) _ f f ^, 1 ^5G'- 1 (45) 



sum*) 



+ II G{u) %M G{52)d4d5 (68) 



we obtain 

r(12; 3) = 5(12)5(13) -JJff ^|^ g(56)r(67; 3)G(74)d4 d5 d6 d7 (69) 
and analogously 

P(12) = -ifi^^l = - ih H G ( 13 ) r (34; 2)G(41)c(3 c(4 . (70) 

The Eqs. (58), (66), (61), (69) and (70) constitute Hedin's set of integro-differential equa- 
tions, whose self-consistent solution, in principle, solves the many-electron problem ex- 
actly. Unfortunately, they are not just numerical relations but contain a functional deriva- 
tive in Eq. (69). Therefore, the Hedin equations cannot be solved self-consistently by com- 
puter codes, but they may be iterated analytically in order to derive useful approximations. 
In practice we can only perform one iteration. We start with the Green function Go of the 
noninteracting system, which corresponds to the single-particle Hamiltonian (Eq. (64)). 
As the corresponding self-energy vanishes in this case, the set of equations simplifies to 

T(12; 3) = 5(12)5(13), (71) 
P(12) = -ifiG (12)Go(21), (72) 

W(12) = v(12) + J J v(13)P(34)W(42)d3d4, (73) 

£(12) = t«? (12)W(l + 2) , (74) 

G(12) = G (12) + J J G (13)S(34)G(42)d3d4. (75) 

The polarization function here corresponds to the bubble diagram of Feynman's diagram- 
matic approach to many-body perturbation theory and leads to the random-phase approx- 
imation for the screened interaction (cf. Figure 5). The expression for the self-energy in 
this first iteration coined the name GW approximation. 



Appendix B 
Quasiparticle equation 

Inserting Eq. (12) into the equation of motion for the Green function of a stationary system 
in the frequency domain 

[fiw-fto(r)] G(ry,w)- f E(r, r"; w)G(r", r'; oj)d 3 r" = 5 (r - r') , (76) 
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which is equivalent to the Dyson equation, yields 

^ to- * I in {[^- ho H ^ (r) - / s (r ' r " ; w) ^ (r " } dV ' } = 6 (r - r ' } ' 

(77) 

Now we multiply with (hui — ef) and take the limit ui — > on both sides. If we assume 
that the system is nondegenerate (i.e., all are different), the left-hand side becomes 

lim (foj - ej) V - ^ (r/) . ( [foj - fto (r)l W - / S (r,r"; W )^ (r")dV'l 

= V* (r') | [ej - ho (r)] ^ (r) - | S (r, r"; ej/h) ^ (r") dV' J , (78) 

and the right-hand side becomes 

lim (fed - Cj ) <S (r - r') = . (79) 

Since V>j (r') does not vanish for all r', the expression in the curly brackets must be zero. 
This leads directly to the quasiparticle equation 

h (r) ^ W + j S (r, r"; ej/h) ^ (r") d 3 r" = e^j (r) . (80) 

It remains valid in the degenerate case, which is seen as follows. We assume that 
the solution of Eq. (80) leads to degenerate amplitudes tpj (r) and energies tj . Then we 
introduce an arbitrary perturbation <f>, e.g., an additional external potential in ho(r), that 
breaks the symmetry in such a way that the degeneracy is lifted. For this (nondegenerate) 
case the above proof applies. The validity of Eq. (80) for the degenerate case is then 
established by taking the limit cj) — ► 0. 

References 

1. P. Hohenberg and W. Kohn, Inhomogeneous Electron Gas, Phys. Rev. 136, B864 
(1964). 

2. A. L. Fetter and J. D. Walecka, Quantum Theory of Many -Particle Systems (Dover, 
New York, 2003). 

3. G. D. Mahan, Many-Particle Physics (Plenum, New York, 1990). 

4. L. Hedin, New Method for Calculating the One-Particle Green's Function with Ap- 
plication to the Electron-Gas Problem, Phys. Rev. 139, A796 (1965). 

5. W. Kohn and L. J. Sham, Self-Consistent Equations Including Exchange and Corre- 
lation Effects, Phys. Rev. 140, Al 133 (1965). 

6. L. J. Sham and W. Kohn, One-Particle Properties of an Inhomogeneous Interacting 
Electron Gas, Phys. Rev. 145, 561 (1966). 

7. J. P. Perdew, R. G. Parr, M. Levy, and J. L. Balduz, Density-Functional Theory for 
Fractional Particle Number: Derivative Discontinuities of the Energy, Phys. Rev. 
Lett. 49, 1691 (1982). 



20 



8. W. G. Aulbur, L. Jonsson, and J. W. Wilkins, Quasiparticle calculations in solids in 
Solid State Physics, edited by H. Ehrenreich and F. Spaepen (Academic, New York, 
2000), Vol. 54, p. 1. 

9. P. Duffy, D. P. Chong, M. E. Casida, and D. R. Salahub, Assessment of Kohn- 
Sham density-functional orbitals as approximate Dyson orbitals for the calculation 
of electron-momentum-spectroscopy scattering cross sections, Phys. Rev. A 50, 4707 
(1994). 

10. W. von der Linden and P. Horsch, Precise quasiparticle energies and Hartree-Fock 
bands of semiconductors and insulators, Phys. Rev. B 37, 8351 (1988). 

11. G. E. Engel and B. Farid, Generalized plasmon-pole model and plasmon band struc- 
tures of crystals, Phys. Rev. B 47, 15931 (1993). 

12. M. S. Hybertsen and S. G. Louie, First-Principles Theory of Quasiparticles: Calcu- 
lation of Band Gaps in Semiconductors and Insulators, Phys. Rev. Lett. 55, 1418 
(1985). 

13. R. W. Godby, M. Schliiter, and L. J. Sham, Accurate Exchange-Correlation Potential 
for Silicon and Its Discontinuity on Addition of an Electron, Phys. Rev. Lett. 56, 
2415 (1986). 

14. M. S. Hybertsen and S. G. Louie, Electron correlation in semiconductors and insula- 
tors: Band gaps and quasiparticle energies, Phys. Rev. B 34, 5390 (1986). 

15. R. W. Godby, M. Schliiter, and L. J. Sham, Quasiparticle energies in GaAs andAlAs, 
Phys. Rev. B 35, 4170 (1987). 

16. J. E. Northrup, M. S. Hybertsen, and S. G. Louie, Theory of quasiparticle energies in 
alkali metals, Phys. Rev. Lett. 59, 819 (1987). 

17. R. Keyling, W.-D. Schone, and W. Ekardt, Comparison of the lifetime of excited elec- 
trons in noble metals, Phys. Rev. B 61, 1670 (2000). 

18. G. Onida, L. Reining, and A. Rubio, Electronic excitations: density-functional versus 
many-body Green' s-function approaches, Rev. Mod. Phys. 74, 601 (2002). 

19. M. Petersilka, U. J. Gossmann, and E. K. U. Gross, Excitation Energies from Time- 
Dependent Density-Functional Theory, Phys. Rev. Lett. 76, 1212 (1996). 

20. L. Hedin and S. Lundqvist, Effects of Electron-Electron and Electron-Phonon Inter- 
actions on the One-Electron States of Solids in Solid State Physics, edited by F. Seitz, 
D. Turnbull and H. Ehrenreich (Academic, New York, 1969), Vol. 23, p. 1. 



21 



